Electromagnetic Response of a Vortex in Layered Superconductors 
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We calculate the response of a vortex core in a layered superconductor to ac electromagnetic 
fields with frequencies lj < 2A/h. In this frequency range the response is dominated by order para- 
meter collective modes which are coupled to the vortex-core bound states of Caroli, de Gennes and 
Matricon. Our calculations show that a vortex core has a more complex and richer dynamics than 
predicted by previous theories. The ac field drives an oscillating, nearly homogeneous supercurrent 
in the direction of the electric field, superimposed with a dissipative current flow which has a 
dipolar spatial structure. The order parameter response at low frequencies is an approximately rigid 
collective motion of the vortex structure perpendicular to the external field. This structure becomes 
strongly deformed at frequencies of order i^>0.5A/?i. Coupling of the vortex-core bound states to 
collective modes of the order parameter at low frequencies leads to substantial enhancement of the 
dissipation near the vortex-center well above that of the normal-state. 



Vortex cores play a key role in dissipation processes of 
superconductors in the Abrikosov phase. Bardeen and 
Stephentl modeled the vortex core as a region of nor- 
mal metal in which the excitations in the core respond 
to an electromagnetic field like electrons in the normal 
metallic state. This is a good approximation for dirty 
superconductors with a mean free path, £, much smaller 
than the coherence length ^o- However, in clean super- 
conductors the low- lying excitations in the corCj-are the 
bound states of Caroli, de Gennes and Matricon.a These 
excitations have superconducting and normal proper- 
ties. They are the source of_circulating supercurrents in 
the equilibrium vortex core,0 and they are stroagly cou- 
pled to the condensate by Andreev processes.aa Their 
response to an electromagnetic field is radically differ- 
ent from that of normal electrons. For vortex cores 
one has two fundamentally different origins of dissipa- 
tion. One is dissipation by the collective motion of the 
whole condensate, and the second is dissipation by tran- 
sitions between Caroli-de Gennes-Matricon bound states. 
These processes are coupled because of the strong inter- 
action between the condensate and the bound states, and 
requires a self-consistent treatment of condensate and 
bound-state djaiamics. Earlier calculations neglected 
this coupling. BQ We present the first fully self-consistent 
calculation of the current response to an ac electric field, 
with frequencies comparable with the gap frequency A/h, 
in the vortex core. Our results show that self-consistency 
is essential for even a qualitative understanding of low- 
frequency dynamics and dissipation by the core. 

We consider an s-wave superconductor with a ran- 
dom distribution of impurities in a static magnetic field. 
The applied ac electric field, SE{t) = —^dtSA{t), is 
linearly polarized in x-direction, and its wavelength is 
large compared to ^o- We assume an impurity scattering 
rate such that the superconductor is outside the super- 
clean limit; i.e. all bound states are broadened by an 
amount comparable to or larger than the "mini-gap" , 



A'^/EfB We investigate the intermediate clean regime. 
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{Ef/A)^o, where we expect the model of a "nor- 



mal metal core" to break down. In the intermediate clean 
range we can use the quasiclassical theory of Feijmi liquid 
superconductivity- formulated by JililenbergerJ3 Larkin 
and Ovchinnikov,B and Eliashberg,Q which is a powerful 
method for studying non-equilibrium superconductivity. 

The numerical solution of nonequilibrium transport 
problems is greatly simplified by a new formulation of 
the nonequilibrium quasiclassical equations. Two of 
the new transport equations are a generalization [€>f the 
Riccati-type equilibrium equations of Nagato et al. Jl3 and 
Schopohl and Maki.EiJ We present the central equations 
of the theory and refer readers to Ref. |l2| for notation, 
their derivation and the connection to the quasiclassical 
Green's functions. 

The scalar transport equations for the distribution 
functions, 7" '^, 7"'^, a;^, x"^ , which are functions of mo- 
mentum pf, position R, energy e, and time t, are 
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External fields, the superconducting order parameter 
and impurity scattering enter the transport equations as 
driving fields and scattering integrals through the self- 
energies, -yfi^^'^, v^^A,K^ /:^R,A,K^ l^R.A..K_ rpj^g^^ 
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self energies are functionals of the distribution functions 
and must be determined self-consistently. Equations 
(H^), together with the self-energy equations, are nu- 
merically stable and are the basis for efficient numerical 
algorithms for solving nonequilibrium problems in the 
quasiclassical theory of superconductivity. For the ac 
response of the vortex we linearize Eqs. (0^) in the 
external field.E3 The linear response equations require as 
input the self-consistent equilibrium functions, 7^"^ and 

7eq^: the order parameter and the impurity i-matrix for 
a pancake vortex. With these equilibrium solutions we 
then solve the first-order transport equations together 
with the linearized self-energy equations and the charge 
neutrality condition. The linearization is justified as long 
as the relation e|(5E| <c (hio)'^ /A^q holds. 
The results presented below are calculated for a lay- 
ered s-wave superconductor with a cylindrical Fermi 
surface along the c-direction, isotropic Fermi velocity v^: 
in the ab-plane, isotropic pairing interaction and a large 
Ginzburg-Landau parameter, k ^ 1. Impurity scatter- 
ing is taken into account self-consistently in the Born 
approximation. We consider a moderately clean super- 
conductor with a long mean free path, £ = 10S,o, where 
^0 = hvf /2TrkBTc is the coherence length, and choose a 
low temperature, T = 0.3Tc. 




FIG. 1. Local density of states at different distances (up 
to 8.8^0, spacing 0.8^o) from the vortex center. The center 
of the vortex corresponds to the bright filled spectrum. 

Fig. H shows the calculated local density of states of 
an equilibrium vortex at various distances from the core 
center. The important feature at the vortex center is the 
zero-energy bound state, which is broadened by impurity 
scattering into a resonance of width F « O.GTc- At finite 
distances from the vortex center the bound states corre- 
sponding to different impact parameters of quasiclassical 
trajectories form a one-dimensional band, which is broad- 
ened by impurity scattering. The bands widen because 
of the coupling to the vortex flow field with increasing 
distance from the center, and develop Van Hove singu- 
larities at the band edges. 



FIG. 2. Distribution of induced current density (left) 
and corresponding local electric field (right) in the vortex 
core at frequency uj = 0.3A/?i. The external electric field 
(5Et^ (t) = (5Eo cos Lot points into x-direction. The top row 
is the dissipative response (~ cosojf), and the lower row is 
the reactive response (~ sinut). Distances from the center 
extend up to 6.3^0 ■ 

Fig. ^ shows typical results for the current and field 
patterns induced by the external electric field, (5E(^(t) = 
(5Eo cos cjt. The linear response decomposes into an in- 
phase (oc cosujt) and out-of-phase (oc sintjt) response. 
The in-phase current response at position R determines 
the local time-averaged energy transfer between the ex- 
ternal field and the electrons, while the out-of-phase cur- 
rent response is non-dissipative. The in-phase and out- 
of-phase currents show characteristically different flow 
patterns. The non-dissipative current is nearly homo- 
geneous (lower left), whereas the pattern of dissipative 
currents has qualitatively the form of a vortex-antivortex 
pair (upper left), which is consistent with a rigid shift in 
y-direction of the equilibrium current density. Addition- 
ally, this pattern is strongly deformed at higher frequen- 
cies (see Fig. 5| below). The right-hand side of Fig. |^ 
shows the tota electric field, which consists of the homo- 
geneous external field and the internal field due to charge 
fluctuations induced by the external field. Note that the 
induced flcld is predominantly out-of-phase and exceeds 
the external field in the center of the vortex core. For 
higher frequencies, 0.5A ^'huj< 2 A, the induced dipolar 
field evolves from an out-of-phase dipole along the di- 
rection of the applied field to an in-phase dipole oppo- 
site to the applied field, thus reducing the external field 
by roughly one half. The induced field vanishes rapidly 
for frequencies above the gap edge. These features of 
the current response are robust and appear in all our 
calculations, independently of the temperature or mean 
free path. The charge density fiuctuations responsible 
for the induced field are of the order e^6/^Q, where 

6 = e \ (5Eo I ^qA/ [TiLo)'^ is a small dimensionless para- 
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meter that measures the amplitude of the perturbation. 
Thus, only a fraction of an elementary charge accu- 
mulates periodically in time over an area of ^g, which is 
consistent with the condition of local charge neutrality.E^I 
Note that the induced charge in the vortex core resulting 
from particle-hole asymmetry, as discussed recently,L3 is 
of order e(-^)^. The in-plane current and field response 
(upper row of Fig. |^) shows that the local dissipation 
can be either positive ( "hot spots" in which the in-phase 
current is parallel to the electric field), as well as neg- 
ative ("cold spots" in which the dissipative current is 
anti-parallel to the electric field), and is dominated by 
dissipation in the vortex center. The net dissipation is 
positive, and is obtained by integrating j -E over the vor- 
tex. At higher frequencies the dissipation is dominated 
by higher energy bound states, and is less localized in the 
vicinity of the vortex center. 

The electromagnetic response of the vortex-core is due 
to an interplay between collective dynamics of the order 
parameter and the dynamics of the Caroli-de Gennes- 
Matricon bound states. At frequencies below 0.5A/7i 
the order parameter performs a nearly homogeneous os- 
cillation perpendicular to the applied field, with a ve- 
locity that is 90° out of phase with the applied field. 
To show the deviations of the self-consistent order pa- 
rameter from a rigidly oscillating vortex structure, on 
which all prior vortex dynamics calculations are based, 
we introduce the displacement vector (5Ro(t) defined by 
(5A(R,t) = 5Ro(R,i)- VAo(R). The velocity field of the 
order parameter is then (5v(R, i) = 5f5Ro(R, t). 
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FIG. 3. Velocity field of the order parameter for 
u — 1.2A/h, showing strong order parameter deformation. 
The left panel shows the in-phase motion of the velocity 
field, and the right panel shows the out-of-phase velocity. 
Distances from the center extend up to 4.7^0- 

In Fig. ^ we show that the order parameter velocity 
field for uj = 1.2A/?i is strongly deformed. The vortex 
core is not rigid, but oscillates in-phase and perpendicu- 
lar to the applied field. The velocity of the vortex center 
is of the order oivfS, so the response is far from that of a 
stationary vortex. The amplitude of the vortex core oscil- 
lation increases with decreasing u. However, for huj > A, 
both the phase and amplitude of the vortex core oscilla- 
tion decrease and approach zero above hcu ^ 2 A. 

The coupling between order parameter response and 
the bound states is demonstrated in Fig. ^, which com- 
pares the nonequilibrium spectral current density (re- 



tarded Green's function, averaged with the weight u/x 
over the Fermi surface) obtained from a self-consistent 
calculation with that obtained from a non-self-consistent 
calculation that takes into account the field-induced tran- 
sitions within the band of Caroli-de Gennes-Matricon 
bound states, but freezes the order parameter degrees 
of freedom. The corresponding contribution to the cur- 
rent density is obtained by multiplying these functions 
by — ^tanh^^^/2 and integrating over e.lll 
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FIG. 4. Spectra for the dissipative part of the current 
density, {vfxSg^), at the vortex center (left) and at 0.8^o 
in x-direction (right). The vertical lines correspond to 
e = uj/2. The solid curves are the fully self-consistent 
spectral response. The dashed lines assume a frozen order 
parameter. 

The dissipative part of the response function shows 
peaks in the spectral current density associated with the 
band of bound states in the vortex core. The spectral 
weight near e = uj/2 does not contribute significantly to 
the current response. Fig. || shows that the bound state 
band is shifted to lower energy by the self-consistent re- 
sponse of the order parameter, leading to an enhanced 
dissipative current. 



FIG. 5. Distribution of induced dissipative current 
density in the vortex core at frequencies uj = l.lA/h (left) 
and 1.5A/h (right). The current density at lj = 1.5A/?i 
is scaled by 2.5 relative to the current density at l.lA/?i. 
Higher energy transitions between bound states produce 
dissipative currents at larger distances from the vortex cen- 
ter. Distances from the center extend up to 6.3^0- 

The dissipative currents show nontrivial structure for 
frequencies in the range A < ?ia; < 2 A. Fig. || shows that 
the dissipative currents flow predominantly in regions 
where bound states with energy e^^ — hio/2 are local- 
ized. We interpret these structures in terms of impurity- 
mediated transitions between the Van Hove band edges 
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(shown in Fig. of the bound-state bands of the equiHb- 
rium vortex. The transition rate increases with decreas- 
ing frequency, and for hui comparable with the width 
of the zero energy bound state the dissipation is deter- 
mined by the spectral dynamics of the zero energy bound 
states at the vortex center. The spectral response of these 
states (Fig. ||) leads to a dramatically enhanced dissipa- 
tion near the vortex center, which is much larger than 
the normal-state dissipation. 
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FIG. 6. Local conductivity of a vortex as a function of 
frequency. The lower set of curves shows uj Kea^x as a func- 
tion of distance from the center (shown in the legend) along 
the y direction. The upper set of curves shows ujlmaxx- 
The thick red curves give the real and imaginary parts of 
the conductivity for a homogeneous superconductor with 
the same mean-free-path. The inset shows a comparison 
of the Drude absorption (black dotted) of a normal metal 
with the enhanced absorption at the center of the vortex 
core. Note that ctq = 2e'^Njv] and I = lO^o- 

Fig. ^ shows the conductivity in the vicinity of the 
vortex core as a function of frequency and distance from 
the vortex center. At the vortex center the real part of 
(absorption) increases dramatically at low frequen- 
cies. Further away from the center the absorption is a 
maximum at a frequency corresponding to transitions 
between states near the Van Hove peaks of the bound 
state band. For comparison we show the conductivity 
of a homogeneous superconductor with the same mean- 
free path, which shows an absorption edge at hu> = 2A. 
Also note that the low-frequency absorption in the vortex 
center is much larger than the Drude conductivity of the 
normal state (inset of Fig. |^). The conductivity sum-rule 
is obeyed; the apparent excess weight of J 'Rc(Jxx{'^)d^ 
is compensated by a negative delta function contribu- 
tion associated with counterflowing supercurrents near 
the vortex center. Fig. ^ also shows the enhanced su- 
percurrents in the vortex center, compared to those of a 
homogeneous superconductor, for frequencies compara- 
ble with the width of the states near zero-energy. 



Finally we note that the order parameter response at 
low frequencies is mostly transverse to the electric field 
and that the dissipative current in the core is predom- 
inantly parallel to the applied field., i.e. there are no 
charge currents in y-direction related the order parame- 
ter motion. However, there is a substantial energy flux in 
the vortex core. Our calculations show, that the energy 
flux is predominantly in direction of the order param- 
eter oscillation and confined in the vortex core within 
a few coherence lengths. The dipolar form of the dis- 
sipative current pattern leads to backflow currents far 
apart from the vortex core, which result locally in "cold 
spots". Thus, the states in the vortex core extract en- 
ergy from the external field, and transport this energy 
several coherence lengths away from the vortex center in 
y-direction. The net absorption is determined by inelas- 
tic processes at the scale of the coherence length. 
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